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Narrowband Biphotons: Generation, 
Manipulation, and Applications 


Chih-Sung Chuu and Shengwang Du 


Abstract In this chapter, we review recent advances in generating narrowband 
biphotons with long coherence time using spontaneous parametric interaction in 
monolithic cavity with cluster effect as well as in cold atoms with electromagneti- 
cally induced transparency. Engineering and manipulating the temporal waveforms 
of these long biphotons provide efficient means for controlling light-matter quan¬ 
tum interaction at the single-photon level. We also review recent experiments using 
temporally long biphotons and single photons. 


1 Introduction 

Entangled photon pairs, termed biphotons, have been benchmark tools in the field 
of quantum optics for testing fundamental quantum mechanics as well as for de¬ 
veloping applications in quantum information technology, including realization of 
the Einstein-Podolsky-Rosen paradox QID, test of violation of Bell’s inequality 
13 01 a, quantum cryptography and key distribution quantum teleportation 
1710, quantum computation O, etc. Traditional methods of producing biphotons 
include spontaneous parametric down conversion (SPDC) (TO] [TT] [T2]| and spon¬ 
taneous four-wave mixing (SEWM) (131 [Hi [151 in nonlinear solid-state materials. 
Photons from these nonlinear processes in free space forward-wave configuration 
usually have broad bandwidth (> THz) and short coherence time (< ps) such that 
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their temporal quantum waveform can not be directly resolved by existing commer¬ 
cial single-photon counters (which have a typical resolution of about ns). 

These broadband biphotons are not suitable for recently proposed protocols of 
long-distance quantum communication and quantum network based on photon-atom 
interaction IHKTtI, because an efficient photon-atom quantum interface requires 
single fiying photons having a bandwidth sufficiently narrower than the atomic res¬ 
onance (typically in MHz). To reduce their bandwidth, optical cavity may be used 
for active filtering mmsiEoiED. However, multiple cavity modes are resonated 
simultaneously due to the broad gain linewidth of the forward-wave parametric in¬ 
teraction and the biphotons are generated in multiple longitudinal modes. Additional 
passive filters locked to the desired mode of the resonant cavity is thus necessary to 
obtain single mode output with reduced generation rate and increased complexity. 

Phase-matching condition plays an important role in determining the photon 
bandwidth in these nonlinear process. The loose constraint of phase matching in 
the forward-wave configuration make it difficult to generate ultranarrow-bandwidth 
biphotons. To overcome this problem, one may use the cluster effect in a monolithic 
cavity with double-pass pumping. An alternative solution is to take the backward- 
wave configuration where the produced paired photons propagate in opposing direc¬ 
tions. In this chapter, we review recent advances in generating narrowband bipho¬ 
tons with long coherence time using forward-wave and backward-wave parametric 
interaction. In particular, we describe the monolithic resonant down-conversion with 
cluster effect (221 EH and the resonant SFWM in cold atoms with electromagneti- 
cally induced transparency (EIT)(24||25][26l. 


2 Monolithic Resonant Parametric Down-Conversion with 
Cluster Effect 

A resonant parametric down-converter can generate single-mode biphotons without 
external filtering when its gain linewidth is narrower than the spacing of adjacent 
resonant modes. For a doubly resonator where the signal and idler fields are si¬ 
multaneously resonant, the mode spacing is determined by the cluster spacing AQc 
between the simultaneous resonant signal the idler modes. The cluster spacing can 
be calculated by solving the following equations Ea, 

± 1 = M{co)An^ +N{(0)Anc, 

M{0}) = [L/{2Kc)]{2[n^{o))+ni{(Oi)] + o}sn"{o}) + o}in"{(Oi)}, 

N{(0) = [L/{nc)][ns-ni + (0sn,{(0s)-(0ini{(0i)], (1) 
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where Us and rii are the refractive indices at the signal and idler frequencies, and • 
and ^ are the first and second frequency derivatives, respectively. 


If we assume that the idler mode spacing Ai is slightly less than the signal mode 
spacing As and the group velocity dispersion is negligible, the cluster spacing can 
also be obtained by multiplying the idler mode spacing by the number of idler modes 
between two doubly resonant modes. 


AQr 




( 2 ) 


For small As — Ai, the cluster spacing will be much larger than the signal or idler 
mode spacing. Moreover, since AQc decreases with the length of cavity, a mono¬ 
lithic doubly resonator has the largest possible cluster spacing and can be used to 
reduce the number of resonant modes within the gain curve. 


For parametric down-conversion of the forward wave type, a monolithic dou¬ 
bly resonator by itself is not enough for generating biphotons in a single longi¬ 
tudinal mode. This is because the gain linewidth of the parametric interaction is 
larger than the cluster spacing, which allows biphotons to be generated in multiple 
longitudinal modes. As an example, we consider a 3-cm long monolithic PPKTP 
crystal pumped by a continuous-wave 532 nm laser. At type-II phase matching, the 
gain linewidth Afsp = 0.885c/(|ni^^ — \L) = 93.2 GHz and the cluster spacing 

AQc = cj |L) = 52.6 GHz where are the group indices at the sig¬ 
nal and idler frequencies and L is the crystal or cavity length. Because Afsp>AQc, 
multimode biphotons are generated and external filters are required to eliminate 
biphotons in unwanted modes to obtain single-mode biphotons. This results in re¬ 
duced biphoton generation rate in additional to increased complexity of the biphoton 
source. 


2.1 Single-mode output 


To obtain single-mode biphotons without additional filtering, it is necessary to in¬ 
crease the mode spacing and reduce the gain linewidth. This can be achieved by 
using double-pass pumping, for example, by depositing a high refiection coating 
at the pump frequency on the output face of the crystal. The double-pass pump¬ 
ing effectively doubles the length of the parametric interaction or reduces the gain 
linewidth by half. The broaden gain linewidth Af^p = 0.443c/(|^^^^ 1^) is thus 

smaller than the cluster spacing AQc = Consider the previous 

example with a double-pass pump. The gain linewidth Af^p = 46.6 GHz is now 
narrower than the cluster spacing = 52.6 GHz, so biphotons can be generated 
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in a single longitudinal mode without the need of additional filtering. 

The principle of single-mode operation is illustrated in Fig.[2a), where the solid 
curve is the gain profile and the vertical lines represent the frequencies of the si¬ 
multaneously resonant signal and idler modes. With one resonant mode aligned to 
the center of the gain curve, the adjacent modes are outside the gain region. For 
comparison, the gain curve of a monolithic resonator with single-pass pumping is 
also shown as a dotted curve where more than one resonant signal and idler pairs 
are present within the gain curve. 




Fig. 1 Single-mode operation is shown as compared to multi-mode operation. The solid and dotted 
curves are the gain profiles (spectral power density) of a doubly resonant forward-wave parametric 
down-converters with double-pass and single-pass pumping, respectively. The vertical red lines 
represent the frequencies of the simultaneously resonant signal and idler pairs with one mode taken 
at the center frequency /o of the gain curve. Quasi-phasematching is chosen for (a) degenerate and 
(b) non-degenerate frequencies. 


As another example, we consider the generation of single-mode non-degenerate 
biphotons. We assume a 3-cm long monolithic PPKTP crystal pumped by a double¬ 
pass 525.5 nm laser. Such source may be useful for applications such as quantum 
repeater ||28l. For example, the signal photons at 795 nm can be stored in nearby 
atomic memories while idler photons at 1550 nm can be sent through a fiber to 
interfere on a distant beam splitter for creating entanglement between two remote 
locations. As shown in Fig.j^b), the gain linewidth Afdp = 98.1 GHz is narrower 
than the cluster spacing AQc = 110.7 GHz and the biphotons are generated in a 
single longitudinal mode. 
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2.2 Experimental realization 


Chuu, Yin, and Harris l2^ have demonstrated a miniature ultrabright source of nar- 
rowband biphotons using a monolithic resonant parametric down-converter with¬ 
out external filtering. A schematic of the experimental setup is shown in Fig. 
The biphoton source uses a nonlinear crystal of which the end faces are spherically 
polished and deposited with a high refiection coating at the signal and idler wave¬ 
lengths to form a monolithic cavity. The high refiection coating on one end face is 
also highly reflective at the pump wavelength for double-pass pumping. The com¬ 
bination of monolithic cavity, double-pass pumping, and type-II phase matching 
allows the generation of single-mode biphotons near degeneracy. When the pump 
is increased above threshold to produce parametric oscillation, the output is single¬ 
mode, as seen by a scanning Fabry Perot interferometer. This suggest that only a 
single cluster is resonant. 



Fig. 2 A monolithic doubly resonant parametric down-converter for generating single-mode 
biphotons. The pump is double-passed by the end face of the monolithic crystal, which has a 
high reflection coating at the signal, idler, and pump wavelengths. The lens is used to collimate 
the signal and idler flelds. The PBS (polarizing beam splitter) separates the signal and idler flelds 
before they are detected by the single-photon detection modules (not shown). 


A typical coincidence measurement of the biphoton wavepacket is shown in 
Fig.[^ The measured curve shows two asymmetric exponential decays possibly due 
to the different reflectivity for orthogonal polarizations. The biphoton correlation 
time is found to be = 17 ns and the biphoton bandwidth is A to = 2;r • 8 MHz. 
Correcting for the collection efficiency, the generation rate is R = l.lOx 10^ 
biphotons/(s mW) which gives a spectral brightness of R/Ao) = 1.34 x 10^ bipho- 
tons/(s MHz mW). 
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Fig. 3 Glauber correlation function of the biphotons. The coincidence counts are measured as a 
function of the time delay between the signal and idler photons at a pump power of 700 jUW. The 
function has two exponential decays with decay constants of 13.29 ± 0.14 ns and 11.33 ± 0.12 ns 
for time delay greater than or less than zero. The biphoton correlation time is 17.07 ± 0.13 ns. 


3 Backward-wave Biphoton Generation 

Parametric down-conversion of the backward-wave type provides an alternative way 
to realize an ultrabright biphoton source without external filtering. As compared to 
conventional forward-wave schemes, such as spontaneous parametric down conver¬ 
sion (SPDC) in nonlinear crystals, the backward geometry allows a much tighter 
constraint of phase matching that leads to substantially narrow linewidth. 


3.1 General Formulism: Free space 

The backward-wave biphoton generation in free space is schematically illustrated in 
Fig.|^ where the generated field di propagates along the -fz direction and d 2 along 
the -z direction. For SPDC with nonlinearity, this geometry can be achieved 
by quasiphase matching 1^ to submicron periodicity and driven by a single pump 
laser beam |[30l[3T]|32l. For SFWM driven by x^'^^ nonlinearity, the phase matching 
condition can be satisfied by aligning two coherent driving laser beams 13^ [261 . 
As compared to conventional forward-wave schemes, the backward-wave geometry 
allows a much tighter constraint of phase matching that leads to substantially narrow 
linewidth. 
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a2 — 


a^ = 

g^gi(/cioZ-a)iot) 


Fig. 4 Schematics of backward-wave biphoton generation from a nonlinearly driven medium in 
free space. 


We describe the spontaneously generated paired photons as quantized field oper¬ 
ators 


ai(z,r) 

d2{z,t) = (3) 

where co/o are the central angular frequencies and ki^ = co/o /c are the wave num¬ 
bers in vacuum. bi{z,t) are their slowly varying envelopes and follow the Fourier 
transform 


bi{z,t) =-^= jbi{z,(o)e ^^’d(0. (4) 

The frequency domain field operators are governed by the following Heisenberg- 
Langevin coupled equations El: 


- + a,(o»--— 

d , . .Ako 

j.+e2«o) + ,^ 


bi{z,(o) = Ki{(o)bl{z,-(o) + Fi{z,(o), 
bl{z,-(0) = K2{(0)bi{z,C0) + {z,-co), 


(5) 


where ai{(o) = -i^Xi{oiio + and g 2 {co) = - «) describe the 

linear propagation effects with the linear susceptibilities Xi^ are the nonlin¬ 

ear parametric coupling coefficients, Ako is the phase mismatching in vacuum, and 
Fi{z, Co) are the Langevin noise operators. Neglecting the Langenvin noise operators 
that contribute to uncorrelated noise photons, the above equations reduce to 


d , , ,Ako 

j-^+g2{C0) + l — 


bi{z,co) = Ki{(o)bl{z,-(o), 
bl{z,-(o) = K 2 {( 0 )bi{z,( 0 ), 


( 6 ) 


which are subject to the boundary conditions [bi (0, co) , bl (0, co')] = [S 2 (T, co) , ^2 ^')] 

5{co — co') and (^j(0, co')bi{0, co)) = {bliL, co')b 2 {L, co)) = 0. The general solution 
to Eq. •HI is given as 
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bi{L,(a) = A{(a)bi{Q,co) + B{co)p 2 {L,-(o), 
bl{L,-(a) = C{(o)bi{0,O))+D{a))bl{L,-a)), (7) 


where 


A(a)) = 


^sinh(2L/2) + 2cosh(2L/2) ’ 


B{co) = 


2ki 

q + Qcoih{QL/2y 


C(ffl) = 
D{co) = 


-2K2 

q + Qcoth(QL/2) ’ 

Qyo!i+g2)L/2 

qsinh(QL/2) + Qcosh(QL/2) 


( 8 ) 


Here^(ft)) = ai{(0) — g 2 {( 0 ) — iAko and Q{(0) = 2 /q^{( 0 ) + 4ki{(o)k2{co). The sin¬ 
gle photon rates can he calculated from 

Ri = {bl{L,t)h{L,t)) = ^ j \B{(0)\'^d(0, 

R2 = (^5(0,t)h2(0,t)) = ^j \C{co)\^d(0. (9) 

The two-photon Glauber correlation function can be determined by 

Gni.^) = {b\{L,t + T)bl{Q,t)b2{0d)k{L,t + x)) = \\if{x)\^+RiR2, ( 10 ) 


where the two-photon wave function is 

= (S2(0,G^i(L,r-fT)) = — / Bi(o)D''io))e~^^^d(0. (11) 

271 J 

Alternatively, the two-photon wave function can also be obtained from 

\\f{x) = {bi{L,t + x)b 2 {Q,t)) = A(ffl)C*((o)e“'®Vto. (12) 

When the system is conservative and the commutation relation of the field operators 
is preserved, Eq s. (pT] ) and (T^ are equivalent. When the Langevin noise fluctua¬ 
tions exist, Eqs. (|ll|) and ( p^ ) become approximated solutions Ea. The photon pair 
generation rate can be calculated from 

R = jwy:)?dx. (13) 

The normalized two-photon cross-correlation function is 
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- ^12 (^) 

R\R2 




= 1 


R\R2 


(14) 


It is instructive to see what parameters play important role in determining the 
two-photon wave function under some reasonable approximations. To ensure the 
spontaneous photon generation remains below the threshold and muti-photon-pair 
events are suppressed, we work in the low parametric gain regime, | /Ci (co) K ‘2 (m) | ^ 
\q^\. We further take k*i ^ k *2 ^ ic and assume the linear loss and gain of the medium 
to the generated fields can be neglected for the generated fields {i.e. , Xi = Xi)' Under 
these conditions, the ABCD parameters in Eq.([^ can be written approximately 

A ^i^hL^iAkoL/l^ 

B ~ K:Lsinc(zifcL/2)e'^^^/2, 

C ~ -KLsmc{AkL/2)e'^’^^/^, 

D ~ e-i^hLgiAkoL/2^ 


where Aki = XiO^io/i'^c) and Ak = Aki — Ak 2 ^ Ako. The biphoton wave function 
in Eq. {nj then becomes 


Va(t) = ^ [ K(co)^(co)e 
In J 

where 0(m) is the the longitudinal detuning function 


<P{co) = sine 


Ak{co)L 


i[Aki{(0)^Ak2i(0)]L/2 


(16) 


(17) 


Equations and ini' are the same as those obtained following a perturbation 
treatment in the interaction picture 1^ . The biphoton joint spectrum is determined 
by two factors: the nonlinear parametric coupling coefficient k{co) and the longi¬ 
tudinal detuning function 0((o) from the linear phase-matching and propagation 
effects. 


3.2 General Formalism: Resonant SPDC 


Doubly resonant optical cavity plays an important role in obtaining single-mode 
output for ultrabright biphoton generation. Here we develop the theory of resonant 
biphoton generation in the Heisenberg picture. 

We denote the signal and idler fields internal to the cavity by the standing-wave 
cavity operators 


as{t,z) = bs{t)cxp{-i£2ql)sm{kqz) 
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(18) 


where bs{t) and bi{t) are the slowly varying envelopes, Q.q and Qy are the cold 
cavity frequencies, kq = q^jL and ky = rn/L. We assume that only one pair 
of signal and idler fields is resonant with the ^th and rth cavity modes simul¬ 
taneously (we will justify this in the following section). Since only the compo¬ 
nents nonorthogonal to the ^th and rth cavity mode interact with the signal and 
idler fields respectively, the generated dipole moment operators for the signal 
and idler fields, which are proportional to bi{t) exp(—exp(/^^z) exp{iAkz) and 
bs{t) exp{—iQyt) exp{ikyz) exp{iAkz) where the phase mismatch Ak = kp — ky — kq, 
are projected against sin(^^z) and sin(^^z). 

Using the input-output coupling formalism (361 [37l, the equations for the evolu¬ 
tion of bs{t) and bi{t) and their relation to the incident fields are 



(19) 


where bf{t) and bf^{t) are the fields incident on the resonant cavity and K is the cou¬ 
pling constant. With {As^Af} and {r 5 ,r/} denoting the spacing of the cavity modes 
and the mirror refiectivity, respectively, the output coupling rates of the signal and 
idler fields are 7 ^ = zl^(l — r^) and Ji = At{1 — n). With and ^ defined as the 
single-pass power loss for the signal and idler fields in the crystal, the total cavity 
decay rates are + Ys and J] = 2^iAi -f 7 . 

The slowly varying output fields b^^^{t) and are related to the internal 

and incident fields by 


br{t) = vfsbs{t)-bf{t) 
bT\t) = ^ib]{t)-i}^\t). 


( 20 ) 


They can be solved by transforming the coupled equations to the frequency domain 
with the Fourier pair 


b(t) = / Z?((jL)')exp(—/(jL)'^)J(jL)' 


b{t')Qx^{im')dt' 

-00 

and converting to rapidly varying quantities 





( 21 ) 


{^s,i ^q,r) 


(22) 
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The output fields a^^^{co) and a^^^^{—COi) can then be expressed in terms of incident 
fields af{(o) and 


= A(a)) a'^{(0)+B{(0) 

= C{(0) +D{(0) (23) 

with commutators [a^f{( 0 ]),df{( 02 )\ = = ^5jkd{(Oi -CO 2 ). 

For small gain the coefficients are given by 


A(m) 

5(m) 

C(m) 

Z)(m) 


rs/2-i{(0-£lq) 

_ -iK^/Wi _ 

[I]/2 - i{(0 - Qq)] [17/2 + i{(Oi - Q.r)] 

__ 

[I]/2 - i{(0 - Qq)] [17/2 + i{(Oi - Q.r)] 

ri/2 + i((0i- £2r) 


(24) 


With Eq. 
quency 


and (24), we obtain the spectral power density at the signal fre- 


S{co) 




__ 

n[A {(0 - Qqf +r/] [4{(0i - Qr)^ + 172 ] ’ 


(25) 


with the biphoton linewidth Ao) = [( —T}^)/2]^^^, and the 
total paired count rate at exact phase matching, 



r,J7(r,+i7)' 


(26) 


Compared to a non-resonant forward-wave SPDC of the same crystal length and 
pumping power [Eq. (p3))], the generation rate of a lossless resonant backward-wave 
SPDC is thus increased by a factor of 


8^ |v^-v/| 

;rrV2 (v^ + v/) ’ 


(27) 


where ^ is the cavity finesse, r is the mirror reflectivity, and Vs^i are the group ve¬ 
locities of the signal and idler photons. The spectral brightness is increased by a 
factor of rih = rirAcOcIAco with AcOq being the gain linewidth of the parametric 
interaction. 
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If the generation rate of biphotons is small as compared to the inverse of the 
temporal length of the biphoton, the accidental two-photon events may be neglected 
and the time domain Glauber correlation function (namely the biphoton wavepacket) 
is given by 




(28) 


where t = ti — ts is the time delay between the arrival of the signal and idler photons. 
The biphoton correlation time is then Tc = (ln2)(l/f^ -f l/J]). The asymmetry of 
biphoton wavepacket in T is due to the order of detection of the signal and idler 
photons. 


3.3 Single-mode output 

The narrow gain linewidth of the backward-wave parametric interaction is advanta¬ 
geous for generating single-mode biphotons in a monolithic doubly resonant down- 
converter. Consider a 3-cm long PPKTP crystal placed inside a cavity of the same 
length. The cluster spacing AQc\ = In • 1.75 cm“^ and the gain linewidth of the 
parametric interaction A cOg = 271 • 0.08 cm“^ Since AQ.c\ is much broader than 
AcDg, there will only be one doubly resonant mode within the gain linewidth when 
the cavity is properly tuned. This is shown in Fig. where the gain curves of the 
forward-wave type (red) is also shown for comparison. 

The properties of the single-mode biphotons can be calculated as for the case 
of forward-wave interaction. We assume a cavity finesse of 1000. The time-domain 
biphoton wavepacket is plotted in Fig.j^a), which is slightly asymmetric due to the 
different ring-down times at the signal and idler frequencies. The biphoton correla¬ 
tion time is 68 ns. The spectral power density of the biphotons is plotted in Fig.j^b) 
and has a Lorentzian shape. The biphoton linewidth is 2;r • 2.1 MHz. For exact phase 
matching, the spectral brightness of the generated biphotons is 8.16 • 10^ s“^MHz“^ 
per mW of pump power and, as compared to a non-resonant forward-wave source 
of the same material, pumping power, and length, is about 80 000 times higher. 


3.4 Experimental challenge 


The backward-wave parametric interaction was proposed in the 1960’s for imple¬ 
menting a mirrorless oscillator in the infrared regime (391. The experimental demon¬ 
stration, however, was not realized until forty years later f^2\ because of the required 
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Fig. 5 Gain profile of backward-wave parametric interaction (solid) as compared to that of a 
forward-wave interaction (dash). The vertical line at 0 denotes a mode pair of the resonant cavity 
that is taken at the degenerate frequency. The next mode pair, separated by the cluster spacing is 
outside the gain linewidth of the backward-wave parametric interaction. 




Fig. 6 (a) Time domain biphoton wavepacket of the backward-wave spontaneous down-conversion 
within a resonant cavity, (b) Spectral power density at the signal frequency. 


quasi-phase matching with sub-micron periodicity 1301 [311 . 

To construct a backward-wave biphoton source as described above, a 532 nm 
laser may be used as the pump source to generate biphotons at the degenerate 
frequency of 1.064 jum. For type-II phase matching, a KTP crystal periodically 
poled with a periodicity of A = 872 nm is required to accomplish the third-order 
(m = 3) quasi phasematching so that kp = Kq + — h, where the lattice k-vector 

Kg = 2nmlA. Realization of such source thus needs a KTP crystal that is periodi- 
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cally poled with a sub-micron periodicity. Although challenging, it can be done with 
current structuring technique . 


4 Spontaneous Four-Wave Mixing with Electromagnetically 
Induced Transparency 



Fig. 7 Biphoton generation via SFWM in a four-level double-A atomic system, (a) The atomic 
energy level diagram, (b) The backward-wave biphoton generation geometry. 


Here we describe narrowband biphoton generation via SFWM with FIT ||40l |4T1 
in cold atoms. The modeled four-level double-A atomic system and backward-wave 
configuration are shown in Fig.|^ 11) and |2) are two long-lived ground states (such 
as the two hyperfine ground states of alkali atoms) and between them there is no 
electric dipole transition. |3) and |4) are two excited states (in some scheme they 
can be the same state). A pump laser (cOp) excites atoms at the transition |l)^|4) 
with a detuning Ap. A coupling laser (cOc) dresses the states |2) and |3) resonantly. 
In the presence of continuous-wave counter-propagating pump and coupling lasers, 
phase-matched backward paired Stokes (U)^) and anti-Stokes (co^^) photons are spon¬ 
taneously produced following the transitions |4) ^ |2) and |3) ^ |1), respectively. 
We assume the pump laser is far detuned from the transition |1) -> |4) and its ex¬ 
citation is weak such that the majority of atomic population remains in the ground 
state |1). Under this ground-state approximation, we further assume that both the 
pump and coupling laser beams are undepleted in the atomic medium. As shown 
in|7Ja), the coupling laser and the weak generated anti-Stokes field form a standard 
three-level A FIT scheme. As a result, the coupling laser not only participates in 
the SFWM nonlinear process but also renders the medium transparent for the res¬ 
onantly generated anti-Stokes photons. This FIT resonance dramatically enhances 
the SFWM nonlinear photon conversion efficiency. Moreover, in this system, the 
anti-Stokes photons propagate with a slow group velocity due to the FIT effect ll42l 
while the Stokes photons travel nearly with the speed of light in vacuum. Below 
we show that this group velocity mismatching is the key to manipulating the phase¬ 
matching spectrum and thus the biphoton bandwidth. In the backward-wave con- 
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figuration paired Stokes and anti-Stokes photons can propagate collinearly with the 
pump and coupling beams, and can also propagate in a right-angle geometry, de¬ 
pending on how well the relationship = 0 is satisfied, where k^^c are wave 

vectors of the pump and coupling fields. 

The cold atoms (without Doppler broadening effect) are confined within a long, 
thin cylindrical volume of a length L and atomic density is N. The experiments re¬ 
viewed in this chapter are mostly done with cold atoms in a two-dimensional ( 2 D) 
magneto-optical trap (MOT) da. To adapt the theory in Sec. | 3 . 1 [ we replace the 
field index 1 with as and 2 with 5'. The nonlinear parametric coupling coefficients 
{Kas and Ks) for the SFWM process are connected to the third-order nonlinear sus- 
ceptibilities {Xas and 

Kis{(0) = i^^Xa){(Oasa + (0)EpEc, 

Ks{(0) = i^X^^'’*{(OasO - (o)e;e:. ( 29 ) 

Here Ep and Ec are the electric field amplitude of the pump and coupling laser 
beams, resectively. Under the ground-state approximation, the third-order nonlinear 
susceptibility for the generated anti-Stokes and Stokes fields are 


Xas^ {(^asO tt)) 


_ A^AH3M32M24M4i/(£0^^) _ 

{Ap + /714) [\nc p - 4 (ft) + irn){(0 + iyn)] 


4(Ap -b /7 i4)((0 — Elel2-\- iYe){0) + ije) 


( 30 ) 


and 


NHi3H32l^2mi/{£ofi^) _ 

^ ^ (Zl;, + / 7 l 4 )[|^ 2 ,p- 4 (c 0 -/ri 3 )( 03 -/ 7 l 2 )] 

^ _ -AfMl3M32M24M4l/(£0^^) _ 

4 {Ap + i 7 i 4 )(w - ^ 2^/2 - iYe){ 0 ) + ^ 2^/2 - ije)' 


(31) 


jlij are the electric dipole matrix elements, Qc = is the coupling Rabi fre¬ 

quency, and Yij are dephasing rates, respectively. Ap = cOp — CO41 is the pump de¬ 
tuning from the atomic transition | 1 ) ^ | 4 ). CO = (Oas — COasO is the detuning of the 
anti-Stokes photons from the transition | 1 ) ^ | 3 ), and we take COasO = (5)31 as the 




anti-Stokes central frequency. From Eqs. ( | 30 | ) and (| 3 l|), we have Xal\(Oas0 + to) = 

(3) '—' — 

Xs {(^asO + to). As the pump laser is far detuned (Ap 714), we further obtain 


to) ^ *(tOa^o - to). ^2^ = A/|r2cp-(7i3-7i2)^ is the effective cou- 


Xas ^ ( (^asO ' 

pling Rabi frequency. Ye = (7i2 + 7i3)/2 is the effective dephasing rate. The third- 
order nonlinear susceptibility has two resonances separated by Qe and each is as- 
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sociated with a linewidth of 27 ^. These two resonances here indicate two SFWM 
paths. In one path the frequency of the anti-Stokes photons is COasO + ^2^/2 and the 
frequency of the correlated Stokes photons is (O 42 Ap — Qe/'^- In the other path 
the anti-Stokes photons have frequency at COasO — ^2^/2 while the paired Stokes pho¬ 
tons at CO 42 Ap Consequently, the interference between these two types 

of biphotons will appear in the two-photon temporal correlation, as shall be dis¬ 
cussed in Sec. |4.1[ The results obtained here agree with the dressed-state picture 

The linear susceptibilities at the anti-Stokes and Stokes frequencies are, respec¬ 
tively. 






4N\fii3\^{co + iyi2)/{£oH) 

\Qc?--A{(0 + iY\3){(0 + iyny 

jV|At24p(tO-iyi3)/(£o^) 

|r 2 c |2 -4(0) - /7i3)(a) - iyn) Aj + Tm ’ 


(32) 


(33) 


where Qp = fluEpjh is the pump Rabi frequency. The complex wave numbers of 
Stokes and anti-Stokes photons are obtained from the relations ks = {o)s/c)^J\ -^Xs 
and has = {(Oas/c)^l -^Xas^ where the imaginary parts indicate the Raman gain and 
HIT loss, respectively. In an ideal HIT system with zero ground-state dephasing, i.e., 
with 7 i 2 = 0 , the linear susceptibility is Xas{(^aso)=^^ implying zero linear absorption 
of the anti-Stokes photons. This allows the nonlinear optics occurring on atomic 
resonance without absorption and hence enhances the efficiency of the nonlinear 
interaction. 

Taking \Qp\ Ap, Eqs. (32) and (33) give has — kasO + and Xs — ^ 

so that the wave-number mismatching is approximately Ak Co/Vg + ia. Here kasO 
is the central wave number of the anti-Stokes field, Vg is its group velocity, and 
a, the imaginary part of the anti-Stokes wave number, characterizes the HIT finite 
loss caused by the non-zero ground-state dephasing rate 712 . Now Eq. can be 
approximated as 


0(m) ^ sine 


coL 

Wg 





(34) 


In Eq. (34), a = 2A^(Ti37i27i3/(|^cr+ 4712713 ) where (T 13 = 2;r|Aii3 p/(£ 0 ^+ 13713 ) 


is the on-resonance absorption cross section of the transition |1) ^ |3). The longi¬ 
tudinal detuning function of Eq. ( [34| ) has a full-width-at-half-maximum (EWHM) 
phase-matched bandwidth determined by the sine function, AcOg = 2;r x 0.88/Tg, 
where Tg = LjVg is the anti-Stokes group delay time. The group delay time can 
be estimated from Tg = LjVg { 2 yi 3 /\Qc\^)OD with the optical depth defined as 


OD = NG 13 L. 

Therefore, there are two important characteristic frequencies that determine the 
shape of the biphoton waveform. The first is the coupling effective Rabi frequency 
Qe , which determines the two-resonance spectrum of the nonlinear susceptibility. 
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The second is phase-matching bandwidth A cOg. In the time domain, they correspond 
to the Rabi time and the group delay time Tg. The competition between 

Tr and Tg will determine which effect plays a dominant role in governing the feature 
of the two-photon correlation. Therefore, in the following we will discuss the two- 
photon joint-detection measurement in two regimes, damped Rabi oscillation and 
group delay. 



100 


Fig. 8 Two-photon correlation function in the damped Rabi oscillation regime. The missing of 
other oscillation periods in (B) is due to the short dephasing time of about 33 ns. The experimental 
parameters used here are 713 = 714 = 27r x 3 MHz, 712 = O. 6713 , OD =11, Ap = — 7 . 5713 , and 
Qp = 0 . 8713 . Data are taken from Ref. (33]. 


4.1 Damped Rabi Oscillation Regime 

We first look at the regime where the the nonlinear spectrum plays a dominant 
role and the optical properties of the two-photon amplitude ( p^ are mainly deter¬ 
mined by the nonlinear coupling coefficient. This regime requires that the effective 
coupling Rabi frequency Qe be smaller than the phase-matching bandwidth A(Og, 
i.e.,Qe < A COg, or equivalently, Tr > Tg, so that we can treat the longitudinal detun¬ 
ing function as 0((o) 1. The biphoton spectral generation rate is proportional to 

I /cLp oc Hence both biphoton spectrum intensity and emission rate are pro¬ 

portional to OD^. The two-photon temporal correlation function exhibits a damped 
Rabi oscillation resulting from the interference between the two resonances of the 
nonlinear coupling coefficient. 

We consider the case of Q.c > |7i3 ~ 7i2| first, which implies a real effective 
coupling Rabi frequency Following Eq.(p^, the two-photon wave amplitude 
now is determined by the Fourier transform of the nonlinear coupling coefficient 
(29) to give 




( 35 ) 
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and 0 (t) is the Heaviside step function, i.e., 


D _ : Nlli311321^241^41 V^as^s 

were ^ ^ceoh^neiAp^iyu) 

0(t) = 1 for T > 0, and 0 (t) =0 for t < 0. The physics of Eq. (35) is under¬ 
stood as follows. Because the two-photon state is entangled, it cannot be factor¬ 
ized into a function of tas times a function of ji/rp, depends only on the relative 
time delay T, which implies that the pair is randomly generated at any time. The 
first term in the bracket on the RHS of Eq. ( [^ represents the two-photon am¬ 
plitude of paired anti-Stokes at cOasO + *^2^/2 and Stokes at cOsO — *^2^/2; while the 
second term is the two-photon amplitude of paired anti-Stokes at COasO — ^2^/2 and 
Stokes at COso + ^2^/2. Equivalently, this frequency entangled state can be written 
as I COasO + ^el^)as \ COsO “ ^e/^)s “ | COasO “ ^el^)as \ + ^el^)s^ To further See the 
interference, let us look at the two-photon Glauber correlation function. 


g(2)(t) = ;^|BL|V^^^^[l-cos(i2^T)]0(T), 


1 


(36) 


which displays a damped Rabi oscillation with an oscillation period of In jQ^e and a 
damping rate of 2/^. The heaviside function 0 (t) shows that the anti-Stokes photon 
is always generated after its paired Stokes photon by following the EWM path shown 
in Eig.|7Ja). The two-photon correlation function also shows the well-known anti¬ 
bunching effect [G*^^^(0) < Similarly to the polarization entangled Bell 

states, the visibility of the Rabi oscillation, resulting from the two-photon interfer¬ 
ence in time domain, can be taken as an evidence for the time-frequency entangle¬ 
ment. 

Eigure [^a) shows the first experimental demonstration of this type of biphoton 
source by Balic et al. 1^ . Eigure [^a) corresponds to the case of Te > Tf and the 
oscillations are clearly resolved. Eigure [^b) is the case of Tg < Tf where only the 
first oscillation is observable and other disappear because of the fast dephasing rate 

%• 

Rabi oscillations have also been observed in three-level and two-level atomic 
systems El HU 091. 


4.2 Group Delay Regime 

Suggested by Balic et a/. 1331 and demonstrated by Du et a/. 1241 . the group delay 
regime is defined as Tg > Tr and the EIT slow-light effect can be used to control 
the biphoton temporal coherence time. The group delay condition is equivalent to 
Qe > ^COg ^ InjTg, i.e., the biphoton bandwidth is determined by phase matching. 
Eor this reason, in this subsection we treat the third-order nonlinear susceptibility 
as a constant over the phase-matching spectrum. As a consequence, the double res¬ 
onances of biphoton generation are suppressed. The two-photon correlation tends 
to be rectangular shaped, more like that of the conventional SPDC photons. We ob¬ 
tain the biphoton wave function approximately as va(t) ^ KoL<P(t), where Kq is the 
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on-resonance nonlinear coupling constant and = I/{2k) f <P(co)e~^^^dco is 
Fourier transform of the longitudinal detuning function. When the FIT bandwidth 
[AcOtr — |^ 2 cp/( 27 l 3 \/OD)] is wider than the phase-matching bandwidth, the anti- 
Stokes loss is negligible and the two-photon wave function approaches a rectangular 
shape. However, when the FIT loss is significant, the biphoton waveform follows an 
exponential decay. 




Fig. 9 Two-photon correlation function in the group delay regime, (a) Ideal rectangular-shape 
correlation with a relative group delay L/Vg and pair emission rate R. (b) Experimental data (red) 
and theoretical curve (solid blue line) obtained from the Heisenberg-Langevin theory. Data are 
taken from Ref. Ei. 


When AcOtr > AcOg, we ignore the FIT linear loss and rewrite Fq. (34) as 



\COL\ 


'wj 


We then obtain the biphoton wave function 

i/r(T) KoL<P{t) = KoVgn{T;0,L/Vg)e~^^^^^. 


(37) 


(38) 


The rectangular function is defined as JT(t; Ti, T 2 ) = 0 (t — Ti) — 0 (t — T 2 ). Fquation 
( [38] ) shows that the anti-Stokes photon is always delayed with respect to its paired 
Stokes photon because of the slow light effect. The two-photon correlation time is 
determined by the group delay Tg = L/Vg. Using Fq. (13) we get the photon pair 
generation rate 
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R=\ko\\l. (39) 

Thus in the group delay regime the total rate of paired counts scales linearly as 
OD even though the spectral generation rate, KqL<P{(o), scales as OD^. This is 
because the bandwidth reduces linearly with optical depth os. As illustrated in 
Fig. I^a), the rectangular-shape waveform can be understood in the following pic¬ 
ture. When the photon pair is produced from the front surface (z = L), the anti- 
Stokes photon has no delay and both photons arrive at detectors simultaneously. 
When emitted from the back surface (z = 0), the anti-Stokes photon is delayed 
relative to the Stokes photon by T^. Since the photon pair generation probability 
density is evenly distributed in the medium, a rectangular shape shows up in the 
biphoton temporal waveform. Figure |^a) shows an ideal rectangular-shape correla¬ 
tion with a group delay of Tg = L/Vg. For conventional SPDC photons, the sub-ps 
rectangular-shaped biphoton waveform has only been indirectly confirmed in tSOll . 
For the SFWM in cold atoms discussed here, we find that to observe the rectan¬ 
gular shape, the condition AcOtr > AcOg sets a lower bound for the optical depth. 
Using AcOtr — \Qc\'^/{2Yl3^/OD) and AcOg ^ In/Tg 7r\Qc\'^/{yi^OD), we obtain 
OD > Or equivalently, it requires the FIT delay-bandwidth product greater 
than two ED. 



Fig. 10 Biphotons with coherence time up to about 2 jUs. Operating parameters are: OD=130, 
Qc = 271 X 1.11 MHz and Dp = 271 x 1.14 MHz. Data are taken from Ref. 1^ . 


In Fig.|^b) the experimental demonstration of a near-rectangular-shape correla¬ 
tion is obtained at the optical depth around 53 ||24l. The experimental parameters 
used here are 713 = 714 = 2;r x 3 MHz, 712 = O.O 2713 , OD = 53, Ap = 48. 67713 , 
Qp = 1.1 6713 , and Qc = 4 . 2 O 713 , respectively. The BIT transparency width is esti¬ 
mated around 3.63 MHz and the phase-matching spectral width about 2.93 MHz. It 
is found that the exponential-decay behavior at the tail is due to the finite BIT loss, 
which alters the correlation function away from the ideal rectangular shape. 

The sharp peak shown in the leading edge of the two-photon coincidence counts 
in FigJ^b) is the first observed Sommerfeld-Brillouin precursor at the biphoton 
level 1^ . In the stationary-phase approximation, starting from Eq. one can 
show that the sharp peak results from the beating of biphotons which are generated 
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outside of the BIT opacity window. The detailed analysis of the precursor generation 
at the two-photon level has been presented in l(52l . 

Following Eq. (38), one expects the length LjV^ of the rectangular function in¬ 
creases if one reduces the anti-Stokes group velocity Vg by reducing the coupling 
laser power. However, Eq. ( [38] ) is valid only when the FIT is preserved. When the 
EIT loss becomes significant, the biphoton has an exponential decay waveform 1^ . 
Therefore at a certain OD, there is a limitation to prolonging the coherence time 
while maintaining the rectangular-like shape. Most recently, by pushing the OD 
to 130, the two-photon coherence time was extended to nearly 2 /is, as shown in 
Fig. [TO which corresponds to a FWHM bandwidth of 0.43 MHz ESll . 


5 Manipulation of Narrowband Single Photons 


The long coherence time of these narrow-band biphotons allows us to not only di¬ 
rectly resolve their temporal waveform with existing commercial single photon de¬ 
tectors (with a typical resolution of 1 ns) but also manipulate their quantum wave¬ 
form with external phase-amplitude modulators ||53l[54l. In addition, compared with 
SPDC using a single pump laser, there are more freedoms to manipulate the the 
biphoton generation in SFWM with two driving lasers. We will show in this section 
how to engineer the biphoton quantum states by controlling the temporal (SS] [56l 
and spatial patterns Enin as well as the polarizations ESI Ell of the classical 
driving fields. 

1) Electro-Optical Modulation of Heralded Single Photons 



Fig. 11 Schematic of heralded single photon generation and conditional modulation with an EOM. 


It is well know that time-frequency entangled photon pairs can be used to effi¬ 
ciently produce heralded single photons with well defined relative time origin. Here 
we describe how single photons may be modulated so as to produce single-photon 
waveforms whose amplitude and phase are functions of time. As shown in Fig. [TT] 
The detection of a Stokes photon at Di triggers the function generator that drives 
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Vt^ (ns) 

Fig. 12 Waveforms of heralded anti-Stokes photons, (a) Modulated and unmodulated waveforms, 
(b) Waveforms with Gaussian and rising exponential shapes. Data are taken from ED 


the electro-optic modulator (EOM) to shape the waveform of the anti-Stokes pho¬ 
ton. In this way the heralded anti-Stokes photon temporal wave function (phase and 
amplitude) may be shaped in the same manner as one modulates a classical light 
pulse. 

When a Stokes photon is detected at (zs = 0), the heralded single anti-Stokes 
photon wave packet shaped by the EOAM can be described as 


= {0\m{T)das{z,T)ds{0,0)\Ws^as) 

= m{T — z/c)\f/{T — (40) 


where \%^as) is the two-photon time-frequency entangled state 1 ^ . and m(T) is the 
modulation. It is clear that the heralded single-photon waveform is directly modu¬ 
lated by the EOM. Therefore, within the coherence time, the single-photon wave¬ 
form can be arbitrarily shaped. As shown in Eig. 12 The heralded single photons 
can be modulated into two rectangular pulses, Gaussian or time reversed exponential 

31 . 


The importance of the electro-optic method is its speed and ability to modulate 
phase as well as amplitude. The technique provides the technology for studying the 
response of atoms to shaped single photon waveforms on a time-scale comparable 
to the natural linewidth. 















Narrowband Biphotons: Generation, Manipulation, and Applications 


23 


2) Quantum Fourier Cosine Transform: Modulation and Measurement of 
Biphoton Waveform 



Fig. 13 Modulation and measurement of time-energy entangled photons synchronously driven 
sinusoidal modulators and slow detectors. 


As described previously, the Stokes-anti-Stokes two-photon temporal waveform 
is directly measured as coincidence counts between two single photon detectors. 
However, if the speed of the detectors is slow as compared to the two-photon correla¬ 
tion time, the biphoton temporal waveform can not resolved by the resolution of the 
detectors. That is the reason why the wide-band SPDC biphoton waveform can not 
be directly measured by the existing commercial single photon detectors. Here we 
describe an approach to the problem of measurement of biphoton waveforms using 
slow detectors. Figure p3] shows the essential idea. The Stokes and anti-Stokes pho¬ 
tons are incident on synchronously driven sinusoidal amplitude modulators (MODI 
and MOD2). The coincidence count rate between single-photon counting modules 
(SPCMs) is measured as a function of the sinusoidal modulation frequency. The 
SPCMs are slow as compared to the pulse width of the biphoton waveform. With T 
equal to the relative arrival time of the Stokes and anti-Stokes photons, the inverse 
Fourier transform of the (frequency domain) measurement of coincidence counts 
versus frequency yields the Glauber correlation function and therefore the 

square of the absolute value of the biphoton wave function 15411 . 

The modulated correlation function can be written as 

G^^\t,t + z) = |mi(r)p|m2(f + T)pG|,^V)> (41) 

( 2 ) 

where mi and m 2 are the two amplitude modulation functions. Gq (t) is the un- 

( 2 ) 

modulated Glauber correlation function. We average t) over a period T 

of the modulating frequency to form the time averaged correlation function 




( 42 ) 
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Combining Eqs.([4T]) and we obtain 


GmV) 

M(t) = 2 / \mi{t)\^\m2{t + T)\^dt 

1 JO 


(43) 


where M(t) is the intensity correlation function of the modulators. Here we assume 
both channels are modulated by the same sinusoidal amplitude modulation m\{t) = 
m 2 {t) = cos(mr+ (p). Then we haveM(T) = 1/4+1/8cos(2(0 t), and the integrated 
coincidence count becomes 


poo ——- 1 poo 

/ / [2 + cos(2mT)]G(^^(T)JT. 

Jo O Jo 


(44) 


We neglect the dc term and normalize to obtain the Fourier cosine transform pair 


F{2co) 



i: 

poo 

Jo 


cos(2(0t)Jt, 
F {2co) cos{2coT)dco. 


(45) 


Therefore the two-photon correlation function can be obtained by an inverse Fourier 
cosine transform from integrated coincidence counts (by slow detectors). 

The proof of principle of the above Fourier technique in measuring biphoton 
waveform was demonstrated with narrow-band biphotons and the slow detectors 
are simulated with integration El. The result agrees well with the directly mea¬ 
sured waveforms. For practical applications, this technique can be extended to short 
biphotons. 


3) Shaping Biphoton Temporal Waveforms with Temporally Modulated Clas¬ 
sical Fields 




Fig. 14 Schematic of biphoton generation with modulated driving-laser fields, (a) Experimental 
configuration, (b) Atomic energy level diagram for multichannel biphoton generation. 
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By passing long single photons through electro-optical modulators (EOM), it is 
possible to modulate one of the paired photons or their correlation function. How¬ 
ever, these external modulators always introduce losses and attenuations. As shown 
in Fig. [ 7 ] due to the two-photon resonance condition in the SFWM process in the 
double-A atomic system, the central frequency of the Stokes photon follows the 
pump laser, while the central frequency of the anti-Stokes photon follows the cou¬ 
pling laser (in the group delay regime). As a result, if we add a modulation to the 
pump (coupling) laser, the modulation will affect the generated Stokes (anti-Stokes) 
photon. As compared to the SPDC source with a single pump laser, the SFWM 
scheme with two driving lasers allows more freedoms to manipulate the biphoton 
state by controlling the two classical driving lasers. 

When the pump and coupling lasers are periodically modulated. As shown in 
Fig. p^b), the driving laser fields are decomposed in frequency domain into dis¬ 
crete frequency components. As a result, biphoton generation follows many possi¬ 
ble SFWM paths. The interference between these multichannel FWMs provides a 
controllable way to manipulate and engineer the biphoton wave packets. Under the 
condition that the atomic optical depth in the transition |1) |3) is high, the EIT 

slow light effect is significant, and there is no spectral overlap between these SFWM 
channels, following the theory by Du et al. 1551 . the biphoton wave amplitude can 
be obtained analytically 


'i>{ts,tas) = -^^=Ep{ts)Ec{tas)Wo{tas-ts) X £ , 

2jlplc 


(46) 


where Ep{ts) and Ec{tas) are the complex envelopes of the pump and coupling 
laser fields. Ip and Ic are the pump and coupling laser average intensities. — 

is the original two-photon wave packet without modulation on 
the driving fields. As shown in Eq. (|46]l, the two-photon time-frequency entan¬ 
glement information is preserved while the pump and coupling field profiles are 
mapped into the biphoton waveform. The Stokes and anti-Stokes temporal correla¬ 
tion function is 




^p{ts)hitas) ^(2) ^ 

^0 ) 5 




(47) 


where Ip{ts) and Ic{tas) are the pump and coupling laser intensity temporal profiles. 

= |v/b(T)p is the correlation function without modulation. Therefore the 
time-averaged correlation becomes 

R{t)=C{t)Ro{t), (48) 

where C(t) = lim^ t^oo .^7 y ^ h {t)Ic {t + z)dt is the time-averaged pump-coupling 

Zi 1 1 pi ^ 

correlation function. Equations ([46|)-([48]) show that it is possible to manipulate the 
biphoton temporal wave packet and its correlation functions in a controllable way. 
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The proof of principle has been experimentally demonstrated in Ref. (561. 

4) Shaping Biphoton Temporal Waveforms with Spatially Modulated Classical 
Fields 

Following Rubin’s group delay picture ini, in the ideal group delay regime, the 
rectangular waveform reflects the uniform spatial distribution of the photon pair 
generation probability. This is indeed ensured in our previous theoretical model 
where the atomic density, pump and coupling laser fields are assumed to be uni¬ 
form. It is noticed that the biphoton generation probability is proportional to the 
pump field intensity. Therefore, if the pump field is spatially modulated along the 
longitudinal direction of the biphoton generation, this spatial modulation will affect 
the biphoton temporal waveform. This effect was first reported in a recent exper¬ 
iment for achieving the two-photon coherence time up to about 2 /is (25l . In this 
experiment, the transverse Gaussian amplitude profile of the pump laser beam is 
projected to the biphoton longitudinal direction and revealed in the two-photon cor¬ 
relation function. This effect becomes only apparent at a high OD of more than 100 
- at OD<60, the theory in Sec. |3.1| still holds and the pump field can be treated as a 
plane wave even though in reality it is in a gaussian mode. 



Fig. 15 Biphoton temporal waveform reveals the Spatial Gaussian profile of the pump laser beam. 
Operating parameters are: OD=130, = 27r x 11.34 MHz and Qp = 2k x 1.14 MHz. Data are 

taken from Ref. (25] 


In Ref. (251, the pump-coupling laser beams are aligned with am angle of 2.8^ to 
the biphoton longitudinal z-axis and the pump beam transverse Gaussian profile is 
projected to the z-axis. Taking into account the pump field profile effect, we modify 
Eq.([n) to have 

i/a(t) = ^ y , (49) 

where the longitudinal detuning function is replaced by . F{Ak) 

is the Fourier transform of the pump field profile f{z) = 1/(2;r) / dkF{k)e^^^ along 
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the z-axis. In the group delay regime, the spatial phase propagation in Eq. ( [49] ) can 
approximated as 


{kas + ks)L/2 ~ + (OXg/2, (50) 

where 0o is a constant phase factor. The two-photon spectrum is mainly determined 
by the phase-matching longitudinal function F[A^(()L))], and k{(o) ^ Kq varies slowly 
in frequency. The we can reduce Eq.(|4^ to 

xi/{t)c,KoVgf{L/2-Vgr)e‘^. (51) 

It is clear that the pump field spatial variation is mapped into the two-photon quan¬ 
tum temporal waveform with its origin delayed by L/ {2Vg) = Tg/2. The two-photon 
temporal correlation time is determined by the group delay Tg of the slow anti-Stokes 
photon. 

Figure shows the experimental results at OD=130 12511 . There are two main 
features of the two-photon correlation function. The fast oscillating spike at the 
leading edge is the biphoton optical precursor which travels at the speed of light in 
vacuum |[52l . The later, slowly varying long waveform is generated from the narrow 
FIT window. The Gaussian shape reveals the pump laser intensity profile as we 
expected from Eq. 

The detailed discussion of the effects of the nonuniformity of the driving fields 
on the biphoton waveform and biphoton engineering with spatially modulated driv¬ 
ing lasers can be found in |[35]| . 

5) Polarization Entanglement 


The SFWM scheme provides a natural entanglement mechanism in the time- 
frequency domain, but it is extremely difficult to produce polarization entanglement 
because of the polarization selectivity of EIT in a non-polarized atomic medium 
(581. The “writing-reading” technique with optical pumping provides a solution to 
polarization entanglement but results in reducing time-frequency entanglement (5^ . 

Figure [T^shows a schematics for polarization entanglement generation in a right 
angle geometry (481 . A single pump laser is retro-refiected and serve both pump 
and coupling laser beams in the SFWM process. The perfect phase matching con¬ 
dition allows spontaneously generated paired photons to be emitted at right angle. 
To produce entanglement in polarization, we make use of the degenerate Zeeman 
sub states of each hyperfine energy level. By choosing the 2D MOT longitudinal 
symmetry axis as the quantization axis (z-axis), Zeeman states with AMp = 0 are 
coupled by the linearly-polarized pump beams. Conservation of angular moment 
along the z-axis allows two possible circular polarization configurations as shown 
in Fig. 16a): | cr/ G~) and | G~ cr+). From the symmetry of our system, the quantum 
state of the paired Stokes and anti-Stokes photons at the two detectors are described 
by 
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|1>5Si,2(F=2) 


|1>5Si,2(F=2) 



Fig. 16 Polarization entanglement generation in a right-angle SFWM configuration. (a)^^Rb en¬ 
ergy level diagram with two possible polarization configurations for the spontaneously emitted 
photon pairs, (b) Experimental setup with a right-angle geometry. Two sets of quarter-wave plates, 
half-wave plates, and PBSs are inserted for measuring polarization correlation and quantum state 
tomography. The figure is taken from Ref. EH). 



Fig. 17 Experimental setup for producing subnatural-linewidth polarization-entangled photon 
pairs. The polarization entanglement is created by the quantum interference of the two spatially 
symmetric SEWM processes driven by two counter-propagating pump-coupling beams (Li and 
L 2 ). The inserted energy level diagrams are two possible SEWM channels for L\ and L 2 , respec¬ 
tively. The figure is taken from Ref. (57). 


\'l'sAts,tas)) = X ^{\<<yas) + (52) 

However, in this right angle configuration, because of some forbidden transition 
between the Zeeman sub states which can not be coupled by the coupling laser, the 
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anti-Stokes photons do not see a complete EIT. As a result, the system can only 
work at a low OD in the Rabi oscillation regime 1481 . 

Figure shows a more robust scheme for generating polarization entanglement 
for narrowband biphotons in the group delay regime l57l . The pump laser beam is 
equally split into two beams after the first polarization beam splitter (PBSl). These 
two beams, with opposite circular polarizations ((7+ and (7“) after two quarter- 
wave plates, then intersect at the MOT with an angle of ±2.5° to the longitudinal 
axis. Similarly, the two coupling laser beams after PBS2 with opposite circular po¬ 
larizations overlap with the two pump beams from opposite directions. In presence 
of these two pairs of counter-propagating pump-coupling beams, phase-matched 
Stokes and anti-Stokes paired photons are produced along the longitudinal axis. In 
each SFWM path, the polarizations of the Stokes and anti-Stokes photons match 
those of the corresponding pump and coupling fields. 

To obtain the polarization entanglement, we must stabilize the phase difference 
between the two SFWM spatial paths. This is achieved by injecting a reference laser 
beam from the second input of PBS2. The two reference beams split after PBS2 are 
then recombined after PBSl and detected by a photo detector (PD, a half wave 
plate and a PBS are used to obtain the interference), as shown in Fig. [T^ This is 
a standard Mach-Zehnder interferometer to the reference laser. Locking the phase 
difference of the two arms of the Mach-Zehnder interferometer with a feedback 
electronics stabilizes the phase of the two SFWM paths. To avoid its interaction 
with the cold atoms, the reference beams are slightly displaced relative to the pump¬ 
coupling beams but pass through the same optical components. 

By properly choosing the driving laser polarizations and the phase between the 
two SFWM paths, all four polarization-entangled Bell states can be realized for 
subnatural-linewidth biphotons Ell. 


6 Applications 



Fig. 18 Schematics of experimental setup for studying photon-atom quantum interaction. Narrow- 
band paired Stokes and anti-Stokes photons are produced from a cold atomic ensemble in MOTl. 
The anti-Stokes photons pass through an EOM driven by a function generator triggered by the 
detection of Stokes photons at D^. We then send the heralded anti-Stokes photons with amplitude 
modulation to the second cold ^^Rb atomic ensemble at MOT2. The figure is taken from Ref. CD- 
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In the past, experimental manipulation of a single photon interacting with atoms 
have been mostly focused in the frequency domain. Now with the narrowband 
biphoton generation technique described in this chapter, it is possible to control 
the quantum interaction between a single photon and atoms in time domain. In this 
section, we review the applications of heralded single photons whose waveforms 


are shaped by an EOM. A general schematics is illustrated in Fig.[^ We work with 
two-MOT setup. Narrowband Stokes and anti-Stokes biphotons are spontaneously 
produced from cold atoms in the first MOT (MOTl). After detection of a Stokes 
photon, its heralded anti-Stokes photon, after shaping by an EOM, is directed to the 
cold atoms in the second MOT (MOT2). The atoms in MOT2 can be either a three- 
level EIT system with a coupling laser or a two-level system without the coupling 
laser. In the following we show three examples on how the photon-atom interaction 
in MOT2 can be precisely controlled by manipulating the temporal waveform of the 
heralded single photons. 


1) Optical Precursor of a Single Photon 


Now it has been well accepted that both phase velocity and group velocity of 
light in a dispersive medium can exceed c, the speed of light in vacuum ||60l [Mil . 
Then, what is the information velocity of light? Could it be faster than c and vio¬ 
late Einstein’s causality in the special relativity? Answering this question motivated 
early study of optical precursors by Sommerfeld and Brillouin in 1914 ll6^[63l[64ll . 
They showed theoretically that the front of a step-modulated optical pulse (as a car¬ 
rier of 1 bit information) propagating in a dispersive medium always travels at c. 
This front, in the form of a transient wave now known as the Sommerfeld-Brillouin 
precursor, is then followed by the main pulse traveling at its group velocity. Study of 
optical precursors is of great interests not only for fundamental reasons since it is re¬ 
lated to Einstein’s causality, but also for applications because of its connection to the 
maximum speed of optical information transmission ll65l[66l[67l[68l[69l[70ir7n[7^ . 
Now it is clear that the precursor is the fastest part in the propagation of an optical 
pulse even in a superluminal medium |[68l[63l2ll. Classical precursors are entirely 
based on macroscopic electromagnetic wave propagation. 

A single photon is described by quantum mechanics. What is the speed of a sin¬ 
gle photon in a dispersive medium? In quantum mechanics, an observable physical 
quantity, such as the speed, usually takes many possible (discrete or continuous) 
eigenvalues. Some may argue that a single photon event may possibly violate Ein¬ 
stein’s causality while the expectation value of its speed never exceeds c. Although 
The propagation effect of single photons through slow and fast light media has been 
studied previously |7l|7l|75ll26l, all of these experiments focused on the group 
velocity picture. It remains a question whether there is a speed limit for a single 
photon in the quantum nonclassical world. The detection of single-photon precursor 
may shine light to this question as well as the understanding of quantum information 
transmission. 

Optical precursor of a single photon was first observed making use of narrow- 
band heralded single photon with a step-shaped waveform liTTl . with the experi- 
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Fig. 19 Single-photon optical precursors from a step amplitude modulation, (a) The heralded anti- 
Stokes photon waveform with a step modulation, (b) and (c) are two-photon coincidences after the 
anti-Stokes photons passing through the BIT system {Qci = 3.5 713 , OD=20) and two-level system 
{Qc 2 = 0, OD=2.5) in MOT2 respectively. Inset (d) shows a gaussian pulse propagation in the two- 
level system with a peak advancement of about 40 ns (the lower curve) compared to the reference 
pulse (the upper curve). Data are taken from Ref. cn 


mental setup shown in Fig. The modulated heralded anti-Stokes photon wave¬ 
form is shown in Fig. [^a). After the anti-Stokes photon has passed through the 
FIT medium in MOT2 (OD=20), the precursor is clearly seen at the rising edge 
and separated from the delayed main waveform, shown in Fig.p^b). To test single¬ 
photon causal propagation in a superluminal medium, we turn off the coupling laser 
in MOT2 to work in a two-level system, which has a negative group delay, as shown 
in Fig. [igd). we observe a peak advance of at least 40 ns and with about 10% 
transmission compared to the propagation through vacuum. However, as shown in 
Fig.[T^c), there is no observable advancement relative to the rising edge in the sin¬ 
gle photon waveform. That is, using the quantum mechanics interpretation, there 
is no observable probability for single photon traveling faster than c. This result 
indicates that the optical precursor is always the fastest part even in superluminal 
propagation and Einstein’s causality holds. 

The observation of optical precursor of heralded single photons fill confirms 
the speed limit of a single photon in a dispersive medium, which is indeed the speed 
of light in vacuum. The information velocity of a single photon does not follow its 
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group velocity. It also suggests that the causality holds for a single photon. 

2) Optimal Storage and Retrieval of Single-Photon Waveform 

The second application example of narrowband heralded single photons with 
arbitrary waveform shaping is for optical quantum memory. The efficiency of a 
photon-atom quantum interface strongly depends on the temporal shape of single 
photons fTS]! . Although storage of of weak coherent pulses have been demonstrated 
with high efficiencies up to 87% (791 [80l, obtaining such high efficiency with sin¬ 
gle photons remains a technical challenge due to the difficulty in having narrow- 
band single photons with optimal temporal waveform at the operation frequency 

(SBlHllEa. 



Fig. 20 Optimal storage and retrieval of single photons with a storage efficiency of (49±3)%. 
(a) The optimal input (red curve) and output (retrieval, green curve) heralded single-photon wave¬ 
forms. (b) The time-reversed retrieved photon waveform matches the input photon waveform after 
normalization. Data are taken from Ref. (84) 


Making use of heralded narrowband single photons with an optimal waveform 
shaped by EOM, Du et al reported an experimental demonstration of efficient stor¬ 
age and retrieval of narrow-band single-photon waveforms using EIT in a cold 
atomic ensemble (84ll . With the ability to control both single-photon wave pack¬ 
ets and the memory bandwidth, the storage efficiency is obtained up to (49±3)%. 
The experimental result is shown in Eig. To our knowledge, it represents the 
highest storage efficiency for a single-photon waveform to date. Because an effi¬ 
ciency above 50% is necessary to operate a memory for error correction protocols 
in one-way quantum computation (851, this result brings the atomic quantum light- 
matter interface closer to practical quantum information applications. 

3) Coherent Control of Single-Photon Absorption and Reemission in a Two- 
Level Atomic Ensemble 
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When an optical pulse propagates through a dispersive medium, the absorption 
and emission can coherently modify its spectral-temporal components and lead to 
many fundamental and important optical phenomena, such as attenuation, amplifi¬ 
cation, distortion, slow and fast light effects ||40l|42|[86l|87l[88l. For a single photon 
in the absorptive medium, the remission can only occur after the absorption, as re¬ 
quired by the causality. Although this quantum time order has been observed as the 
antibunching effect in resonance fluorescence (491 [89l, the processes in these ex¬ 
periments are chaotic. After a single photon enters a two-level atomic medium and 
gets absorbed, the spontaneously reemitted photon usually goes to 4;r solid angle 
randomly. 

It turns out a directed “spontaneous” emission excited by a single photon is possi¬ 
ble if the timing of absorption is traceable ll90ll . Moreover, this direct “spontaneous” 
emission can be manipulated in time domain by controlling the single-photon tem¬ 
poral waveform ED. For example, making use of the destructive interference be¬ 
tween the emission (or scattering) and the incident photon wave packet, the prob¬ 
ability of remitting the photon during the absorption can be completely suppressed 
when the incident photon has an exponential growth waveform with a time constant 
equal to the excite-state lifetime. The remission process only starts after the incident 
photon waveform is switched off and thus can be controlled on demand. 
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Fig. 21 Coherent control of single-photon absorption and reemission in a two-level laser-cooled 
atomic ensemble. Heralded anti-Stokes photon waveforms after passing through the two-level 
atomic ensemble at (a) OD=0 (vacuum), (b) OD=3, and (c) OD=8. (d) and (e) show the mea- 

( 2 ) 

sured conditional autocorrelation gc for confirmation of the single-photon quantum nature. Data 
are taken from Ref. ED 


We illustrate the process in Fig. 18 where we turn off the coupling laser in MOT2 
and have the atoms there in a two-level system. With the EOM modulation, we 
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produce single anti-Stokes photons with an exponential growth waveform with a 
time constant equal to 1 /( 2713 ) =26.5 ns for the incident anti-Stokes photons, as 
shown in Fig. |^a). Here 2/ = 2;r x 6 MHz is the population decay rate in the 
excite state. Ai t = t 2 — ti =0, the waveform is switched off with a fall time of 
3 ns. After passing through the atoms, the photons are coupled into a SMF to be 
detected by a SPCM. Fig. [^b) shows the result at OD=3. At this modest OD, 
during the exponential growth period (t < 0), the photon waveform is only partially 
absorbed. After the incident photon is switched off at T=0, this partially absorbed 
waveform is released (re-emitted) following a exponential decay curve which is 
determined by the lifetime [1 / (2/)] of the excited state. As we increase the OD, the 


incident photon gets absorbed more heavily. At OD= 8 , as shown in Fig. 21 c), the 
photon is completely absorbed and the probability in finding the remitted photon 
at T < 0 is nearly zero due to the destructive interference. As expected, at T > 0, 
the interference between the incident waveform and the emission disappears and 
we observe the remitted photon. Consequently, the absorption and reemission of the 
single photon is completely separated in time domain. 

This technique can be used to efficiently excite a single quantum absorber in a 
cavity by a single photon ||92j[93l. The result may find potential applications in the 
quantum networks which require efficient conversion between fiying single-photon 
states and local atomic states 1941 . 


4) Single-Photon Differential-Phase-Shift Quantum Key Distribution 


Heralded 



single 



Photons 




{0,n} 


jum 


Alice 



> 


/ 

/ N 

" \ 

\ 

n n 

/ 

^ T 

50:50 50:50 

BsK 


ninin 

c 

SI u 

Dl 

1 

Inin 


Bob 

1 

\ 

\ 

(d) '(c) '(b) 

(a) 


Fig. 22 Schematics of single-photon differential phase shift quantum key generation and distribu¬ 
tion. 


Due to the long coherence time of the heralded single photon, its wave packet can 
be encoded into many time bins with phase-amplitude modulation. It thus become 
an ideal information carrier for the differential phase shift quantum key distribu¬ 
tion (DPS-QKD) ||95]|96l. Traditionally, discrete polarization quantum states have 
been widely used due to their simplicity im. However, the fiber length of such a 
polarization-based QKD system is limited by the birefringence effect that causes 
the polarization fiuctuation on the receiver. This limit can be overcome by the DPS- 
QKD scheme that is polarization independent. The DPS-QKD also shows tolerance 
to photon-number-splitting (PNS) attacks 19611981 . 































Narrowband Biphotons: Generation, Manipulation, and Applications 


35 


A simplified schematics of the first single-photon DPS-QKD system is illustrated 
in Fig. [22| 1991 . Alice divides the single photon into N (> 3) time slots (with a period 
of T). The keys are encoded by the random relative phase shift between consecutive 
pulses in 0 or ;r. Bob detects the incoming photon using an unbalanced M-Z inter¬ 
ferometer setup with a path time delay difference equal to T. Taking A = 3 as an 
example, the detection at Bob’s site yields four possible time-slot outputs, (a)-(d), 
as shown in Fig. As Bob detect a photon, he records the time and which detec¬ 
tor clicks. If the detector clicks at the (b) or (c) time slot. Bob tells Alice only the 
time slot information through a classical channel; otherwise. Bob discards the pho¬ 
ton. Using the time-slot information and her phase encoding records, Alice knows 
which detector clicked at Bob’s site. Defining the clicks at Di and D 2 as “0” and “1” 
respectively, Alice and Bob can obtain a confidential bit string as a sharing key. The 
photon sent from Alice to Bob can be written as one of the following four states: 
(II 1 O 2 O 3 ) ± IO 1 I 2 O 3 ) ± | 0 i 02 l 3 ))/\/ 3 , where li=i, 2,3 represents one photon at time 
slot i. As nonorthogonal with each other, the four states cannot be perfectly iden¬ 
tified by a single measurement, as shown by the noncloning theorem IIOOL which 
guarantees the security of the scheme. 

In the original DPS-QKD proposal (951, a single photon is split into N paths with 
different lengths and then recombined with passive beam splitters. This brings un¬ 
avoidable loss, which results in a low key creation efficiency [oc (A — 1)/A^] when 
A is large. Meanwhile, the phase stabilization between different paths becomes a 
serious problem A increases. In the experimental demonstration of single-photon 
DPS-QKD, each heralded narrowband single photon is modulated into many time 
slots (A up to 15) using a pair of phase and amplitude EOMs (991 . Without usind 
beam splitters, the entire key creation efficiency scales as (A — 1)/A and approaches 
100% at the limit of large A The details of the experiment are described 

in Ref. (93. 


7 Summary 

In summary, we have reviewed recent development in narrowband biphoton gen¬ 
eration. For the SFWM in cold atomic medium, FIT effect is used not only for 
resonantly enhancing the nonlinearity by eliminating the resonance absorption 
but also for tuning the phase-matching bandwidth with its slow-light effect. For the 
monolithic resonant backward-wave and forward-wave SPDC with cluster effect 
and double-pass pumping, it is possible to realize a miniature ultrabright biphoton 
source. These narrowband biphotons with long coherence time from tens nanosec¬ 
onds to several microseconds can be used to produce heralded single photons with 
arbitrarily shaped temporal waveforms by phase-amplitude modulations. We also 
reviewed their applications in manipulating temporal quantum interactions between 
single photons and atoms, as well as in the quantum key distribution. Most recently, 
it was demonstrated that a single photon with time-reversed exponential growth 
waveform can be loaded into a single-sided Fabry Perot cavity with near-unity ef- 
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ficiency I102lll03ll . With the time-resolved quantum-state tomography I1041ll05ll . 
their further applications in quantum network and quantum information processing 
are to be explored. 
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